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ABSTRACT 

In order to understand the mutual influence between active galactic nuclei (AGN) and star 
formation during the evolution of galaxies, we investigate 142 galaxies detected in both X-ray 
and 70|am observations in the COSMOS (Cosmic Evolution Survey) field. All of our data are 
obtained from archive. X-ray point source catalogs from Chandra and XMM-Newton obser- 
vation; far-infrared 70|^m point source catalog from Spitzer-MIPS observation. Although the 
IRAC colors of our samples indicate the existence of star formation, the ratio of rest frame 
2-10 keV luminosity to total infrared luminosity (8-1000 |a.m) shows that AGN predominates 
the spectral energy distribution (SED). We identify obscured AGN in these 70[a.m luminous 
galaxies as characterized by a larger hardness ratio. The higher X-ray obscuration fraction in- 
dicates an extra contribution from the star formation in the host galaxy in addition to the usual 
AGN dusty torus. If AGN feedback occurs in their host galaxies, the star formation must be 
quenched and dust in host galaxy will be dispersed. However, our temperature fitting shows 
there is no significant dust temperature enhancement in far-infrared wavelength, suggesting 
the thermal radiative feedback is absent from the AGN to the galaxy. 

Key words: galaxies: high-redshift - galaxies: active - infrared: galaxies - X-rays: galaxies. 



1 INTRODUCTION 

Both observational and theoretical evidences over the past decades 
strongly suggest the existence of supermassive black hole (SMBH) 
within the spherical bulge of galaxies. The tight connection 
between the SMBH mass and the bulge properties, such as 
bulge luminosity or stellar velocity dispersion, indicates a sce- 
nario of co-evolution fo r centr a l black hole and its host galaxy 
jKormendv & Richston^ 1 19951 : iMagorrian et al.1 [7998). Galaxy 
merger is a plausible interpretation to explain the concomitant 

i growth of black ho le and host galaxy under the ACDM model 
Treister et al. hold) . Such interpretation implies the mergers dis- 
tort galaxy morphology and trigger the intense star f ormation, pro - 
ducing a large amount dust and molecular gas (Magdis et al.l201lh . 
Meanwhile, interstellar medium (ISM) is also funneled onto the 
SMBH in the galactic nucleus and consequently triggers a phase of 
active galactic nuclei (AGNs), with its powerful electromagnetic ra- 
diation spanning from radio to gamma-ray. In the later galaxy evo- 
lution stage, the merger scenario implies that the radiation-pressure 
or kinematic-wind feedback from AGN may quench the star forma- 
tion, disrupting the remnant gas and dust. At the same time, central 
AGN turns into an optically unobscured quasar, a class of luminous 
AGN. Final outcome in galaxy evolution is a red elliptical galaxy 
jHopkins et"ani2008ll . 
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ISanders et al.l ( Il988ah first proposed a connection between the 
Ultra Luminous Infrared Galaxies (ULIRGs; Lm. ^ W^^ Lq) and 
quasars as two successive snapshots of merger event between two 
gas-rich spiral galaxies. Indeed, from IRAS 1 ly ULIRGs stud- 
ies, the merger-ULIRG connection has been r evealed by a large 
fraction of inte racting galaxies at z < 0.1 ( Borne et al .' 2000; 
ICui et alj|200ll ; [Veilleu x et al. 2 002) . Besides, the merger fraction 
has been shown to increase as infrared luminosity raises, merging 
ga laxy bein g prevalent population amongst ULIRGs to redshift ~ 
1 jShi et al . 2009). However at higher redshift ~ 2, the merging 
galaxy fraction drops slightly and non -interacting disk galaxy frac- 
tion increases (Kartaltepe et alj2012l) . 

On the other hand, the merger-AGN connection is more dif- 
ficult to observe. Although, in the local universe, optical luminous 
quasars are considered to be related to post-star burst merger stage 
jCanalizo & Stocktonll200ll ; lBennert et"ai]|2008ll . at 1.5 < z < 2.5, 
AGNs are more likely to occur inside of normal disk or bulge 
galaxy, instead of interacting systems, suggesting that secular pro- 
cesses (such as bars or n ucleus rings) may t rigger nuclear activity 
instead of major-mergers dKocevski et al.l2012) . The potential con- 
nection between ULIRGs and AGNs, is still under investigation, 
and the nature of underlying physical process respo nsible for such 
extreme far-infrared lum inosity is in debate (Sand ers & Mirabell 
ll996l ; lFarrahetal.ll2003h . 

In order to determine the exact role of the AGN in the overall 
galaxy evolution, it is essential to separate the AGN contribution 
from the normal star formation activity. To decouple two compo- 
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nents, we can use three distinct methods: 

(i) Specific emission lines ratio - Fine structure radiation in optical 
and infrared wavelength are reliable tracers of excitation state of 
the ISM in galaxies. A classification based on two independent line 
ratios can help to segregate AGN and star-forming galaxies. For 
instance, investigatin g [OIII]5007A/H/3 versus [NII]6584A/Ha 
( [Baldwin etZlll98ll. BPT) diagram will lead to distinguish the 
photoionization, derived from normal HII regions, from power-law 
continuum generated by AGNs. Unfortunately, galaxies undergo- 
ing merger event suffer from dust obscuration in optical wavelength 
that limits the efficiency of such diagnostics. 

By the advent of the infrared spectrometer (IRS) on the 
Spitzer telescope, additional emission lines were identified for 
AGN diagnostics in the mid-infrared, that are less sensitive to ex- 
tinction. For example, the line ratio [O IV]25.9fJ.m/[Nell]I2.8jim 
or [NeV]I4.3|i.m/[Nell]12.8iim is a reliable indicator to sepa- 
rate AGNs from star-forming galaxies (Annus 2006). The differ- 
ent transition lines of Polycyclic Aromatic Hydrocarbon (PAH) 
molecules are other reliable tracers, wit h weaker PAH emission 
corre sponds to higher AGN contribution (IVoitlll992l : IVeilleux et al.l 
I2OO9I) : 

(ii) Continuum slope - As the influence of the AGN becomes sig- 
nificant in the galaxy, the infrared spectral in dexes start to be domi- 
nated by a power-law contin uum at ~ 10 (im jAlonso-Herrero et al.l 
l200d iDonlev et al.l l2007h . For longer wavelengths, the mid- 
infrared to far-infrared ratio is a simple measure to quantify 
the relative contrib utions of AGNs and star-forming galaxies 
jVeilleux et alj2009l) : 

(iii) Spectral energy distribution fitting - With the development of 
deep multi-wavelength photometric surveys, panchromatic Spec- 
tral Energy Distribution (SED) studies are becoming increasingly 
popular. Theoretical or empirical templates are commonly used 
to fit the SED and determine the photometric redshifts of galax- 
ies. Multi-component fitting, including the contribution from stars, 
dust, and AGNs is now possible thanks to very broad photomet- 
ric coverage. In the case of heavily obscured AGNs, although star 
formation of host galaxy dominates the flux at most wavelengths, 
a bump in mid-infrared produced by the circumnuclear dust radi- 
ation, can be an indicator of an obscured AGN. Therefore, identi- 
fying hot dust component from the galaxy SED is a reliable way 
to determine the AGN contribution to total infrared luminosity 
jMulIanev et alJl201ll ; Fozzi et alj20l'^ . 

Although the different methods may not agree completely on 
defining a pure sample (i.e. pure-AGN or pure starburst galaxy), 
they all provide evidences that, for the major part of the samples, 
AGN and star formation occur concomitantly. Rather than just pur- 
suing the intrinsic AGN contribution fraction on the total infrared 
luminosity, what is essential here is to understand the mutual influ- 
ence on the properties of AGNs and hosts star-forming galaxies. 

The SEDs of Luminous Infrared Galaxies (LIRGs; 10^^ Lq ^ 
LiR < IQ^^Lpi) and ULIRGs are peaked in the range of 40- 
200^m dSanders & Mirabel|[l99 ^. In fact, the 70^m Spitzer band 
is ideal to tmambiguously trace star-forming galaxies as it is lit- 
tle affected by PAH emissions, silicate absorption, and stellar flux. 
Central AGNs are usually identified by their strong X-ray contin- 
uum emission. Indeed, high energy X-ray photons emi tted by hot 
corona of accretion disk (e.g..! Haardt & Maraschill993h around the 
central black-hole are usually little absorbed by dust and gas from 
host galaxies, compared to lower energy UV-photons. One notice- 
able exception are compton thick objects {Nh 1? 10^"* cm~^) for 
which even hard X-ray photons are absorbed. 

Several studies have investigated the nature of the 70 (xm 



galaxy population. IPatel et sA\ ( l201lh lead a spectroscopic follow- 
up in optical wavelength of 70 [xm galaxies selected from the 
Spitzer Wide-area Infrared Extragalactic Legacy Survey. Their re- 
sults suggest that the most of the IR photons are powered by star 
formation, while contribution from AGN dusty torus emission are 
negligible for the non-QSO-dominated samples. From the SED 
study of 61 70^Lm selected gala xies from the . 5 deg ^ wide Ex- 
tended Groth Strip (EGS) field, Svme onidis et al. concluded 
that, even in the presence of powerful hard X-ray emission origi- 
nated from AGNs, dust emission templates are required to explain 
the observed strong far-inf rared luminos i ty. To reveal a potential 
starburst- AGN connection, iTrichas et al.l (|2009|) used a statistical 
K-S test (Kolmogorov-Smimov test) to assess differences in hard- 
ness ratio between X-ray detected 70|J.m galaxies and the whole X- 
ray population. However, shallow observations in both far-infrared 
and X-ray lead to a less significant probability of the K-S test, 
interpreted as low hydrogen column density (Njj) in the nearby 
SMB H. This result is in c ontrast to current galaxy formation sce- 
nario jHopkins et al ]|2008l) . 

In this paper, we aim to dissect the role of AGN in the 
(ultra)infrared luminous phase of galaxy evolution. We extracted 
the X-ray detected sub-sample from unconfuse d 70 [xm catalog. 
Our m ain lQ\im galaxy catalog was published bv lKartaltepe et all 
fcOlOl) . including total infrared luminosity, multi-wavelength pho- 
tometry and redshift. We describe the data used for this work in 
section 2. As we are interested in physical property of AGN and 
far-infrared selected host galaxy, we cross-match X-ray and 70(i.m 
datasets. In section 3, we depict our matching method as well as the 
methods used to estimate and calibrate different physical parameter 
measurements. We present our results in section 4, and then con- 
duct a detailed discussions discussions in section 5. Throughout 
this paper, we adopt a fiducial cosmological model with the fol- 
lowing parameters: Hq = 70 km s^^ Mpc^^, and Qm = 0.3, Aa = 
0.7. Unless otherwise stated, all magnitudes in this paper are in the 
Vega system. 



2 DATA 

The Cosmic Evolution Survey (COSMOS) is the largest treasury 
program using Hubble Space Telescope, imaging a ^2 deg^-wide 
equatorial field in the optical F814W filter (approximately corre- 
sponds to I band) with the Adva nc ed Camera for Surveys (ACS) 
instrument JScoville et al.ll2007allbt iKoekemoer et al.ll2007l) . Sev- 
eral large follow-up campaigns using ground-based and space- 
based telescopes produced the deepest comprehensiv e photomet- 
ric and spectroscopic dataset across a whole spectrum ( Capak et"a l] 



I2OO7I : iHasinger et al.l2007l: ISanders et al .120071: IScott et al.l2008l1H 

Our in itial catalog is a 70\im catalog extracted bv lKartaltepe et al] 
( l2010h . who have publicly released their measured total infrared 
luminosity, redshifts, and multi-wavelength photometry. Here, we 
will give a brief description of our infrared and X-ray datasets. 



2.1 Spitzer-COSMOS 

Spitzer-COSMOS (S-COSMOS) is a Legacy program designed to 
cover the COSMOS field with deep Spitzer observations in IRAC 
four bands (i.e. 3.6, 4.5, 5.8, and 8.0vim) and MIPS three bands 



^ c.f. the COSMOS Special Issue of the Astrophysical Journal Supplement 
Series, in September, 2007. 
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(i.e. 24, 70, and 160(i.m). The initial Cycle 2 program conducted in 
2006 consists of a deep IRAC survey, a shallow MIPS survey, and 
a small deep MIPS survey. The total integration time of the deep 
IRAC survey is 166 hours, reaching 5a sensitivities of 0.9, 1.7, 
11.3, and 14.6 iiiy, in 3.6, 4.5, 5.8, and 8.0(i.m band, respectively 
jSanders et al.l2007h . A complementary MIPS deep observation of 
450 hours has been conducted in Cycle 3, ensuring accurate 70 
and 160iim flux density measurements. The median exposure times 
were ~ 3400s, 1350s, and 270s for 24jim, 70fi.m, and I60(i.m band, 
respectively, corresponding to Scr depths of ~ 0.08, 8.5, and 65 m jy 
jLe Floc'h et alj2009l; iFraver et al.ll2009l; iKartaltepe et"Zll2010l) . 

2.2 X-ray dataset: XMM-Newton & Chandra 

The COSMOS field has been covered with observation from both 
XMM-Newton (XMM) and Chandra telescopes. 

A 2 deg^ contiguous survey has been conducted with XMM 
to reach medium depth (~60 ks). COSMOS X-ray catalog made 
by XMM observation (XMM-COSMOS) provides three different 
energy bands. Base on the logN-logS relationship, the flux limit of 
0.5-2 keV (soft band), 2-10 keV (hard band), and 5-10 keV (full 
band) drop to 7.2 x lO"^'' ergs cm"^ s~\ 4.0 x 10 ~^^ ergs cm~^ 
s~'^, and 9.7 x I0~^^ ergs c m~^ s~^, respectively jHasinger et al.l 
l2007l ; ICapDelluti et alj200 ^. In total, 1887 point-like so urces were 
identified in the XMM-COSMOS survey. iBrusa et al.l ( |2010|) has 
released a cross matched catalog b etween these sour ces and opti- 
cal counterparts (I band, catalog of Capak et al. I (12007,) ). providing 
redshifts (spectroscopic and photometric), UV to 24 [im photome- 
try, and hardness ratio. 

The Chandra-COSMOS (C-COSMOS) Survey images a ~0.9 
deg^ field in the central part of the original COSMOS field, with 
the effective exposure time ~160 ks in center 0.5 deg^ and ~80 
ks in outer 0.4 deg^. C-COSMOS catalog provides three different 
energy bands, 0.5-2 keV (soft band), 2-7 keV (hard band), and 0.5- 
7 keV (full band), with corresponding flux limits of 1.9 x 10^^® 
ergs cm~^ s~^, 7.3 x 10~^^ ergs cm~^ s~^, and 5.7 x 10~^® ergs 
cm~^ s~^, respectively. In total, 1761 point-like source s have been 
extracted in the in C-COSMOS field jSlvis et allioogh . 

In this paper, we combine observations from XMM and Chan- 
dra, taking advantage of the wider observational field of XMM- 
COSMOS and the deeper observational depth of C-COSMOS, to 
maximize the total number of X-ray point source counterparts of 
our IQyaa selected galaxy sample. 



3 METHOD 

3.1 Matched X-ray and far-infrared point source catalog 

The spatial resolution of the Spitzer Telescope at 70(j.m is around 
18", much larger than the spatial resolution of XMM or Chandra 
X-ray observatory. Therefore, i t is crucial to accuratel y detect the 
position of the 70|i.m galaxies. IKartaltepe et alj ( 120101) have iden- 
tified the optical or near-IR counterpart of 70 (im selected galax- 
ies. We took advantage of the Altas resource from IRSA/COSMOS 
archive to match the optically identified counterpart of 70p.m se- 
lected sources with the XMM and Chandra point source catalogs, 
separately. We identified 108 and 92 70 (xm counterparts by match- 
ing with XMM and Chandra using searching radii of 3" and 1", 

^ NASA/IPAC infrared service archive of COSMOS project. 
|http://irsa.ipac.caltech.edu/data/COSMOS| 



respectively. The median angular distances between the position 
of the X-ray source and the position of optical/near-IR identified 
counteipart of our 70iim-selected galaxies are 0.9" and 0.3" for 
the XMM and Chandra point source catalogs, respectively. Fifty- 
eight sources were identified in both X-ray catalogs, and a total 
number of one hundred and forty-two galaxies have been identi- 
fied to have both lQ\im and X-ray measurements in the COSMOS 
field. We will refer to this sample as "70iJ.m/X-ray " galaxy cata- 
log in the remaining p art of this article. From the same datasets, 
IKartaltepe et alj ( I2OI0I) identified 154 X-ray detected 70 ^.m galax- 
ies. The 12 extra objects compared to our own matching catalog 
are originated from a slightly less strict matching criterion. These 
additional sources are more distant from the 70(j.m source position. 
The number of sources from our different catalogs and matched 
catalogs are summarized in Table |T] 

3.2 Redshift determination 

Several spectroscopic surveys were conducted to measure the red- 
shifts of galaxies in the COSMOS field. The most extensive s ur- 
vey on this field is the zCOSMOS survey l lLillvetalJl2007l) , a 
deep spectroscopic survey conducted with the VLT-VIMOS multi- 
object spectrograph to study the evolution of faint/obscured galax- 
ies through direct analysis of their emission lines. Although the 
observations enable us to probe properties of faint objects, only 
the central I deg^ is covered by zCOSMOS. Forty-nine of our 
70iJ.m/X-ray galaxies had their redshifts confirmed spectroscopi- 
cally by zCOSMOS. Observations using other facilities have been 
conducted on this field providing precise redshifts for fifty -three ad- 
ditional 70n.m/X-ray galaxies CAbazaiian et al. .20091 iTrump et al.l 
I2OO7I ; IKartaltepe et alfcoTolJPrescott et alj|2006l) . In total 102/142 
(71.8%) galaxies with X-ray and 70vim detection have spectro- 
scopic redshift. Tabled summarize the origin and number of our 
targets with spectroscopic redshifts. 

For the remaining galaxies, we used the photometric red- 
shift cat alogs built for both non-AGN sources a nd XMM-selected 
sources ('ll bert et al.ll2009t 'Salvato et al."2009'), using the public 
available software Le PHAR^ ( Amoutsetal. 1999; Ilbe rt et all 
[2OO6I) . The 70ixm/X-ray galaxies which do not have spectroscopic 
redshift have been assigned a photometric redshift (40 out of 142 - 
28.2%). 

The photometric and spectroscopic redshifts of 70nm/X-ray 
galaxies follow the overall lQ\im gal axy relation with a scatter of 
0.02 X (1 + 2) ( IKartaltepe et alj20ICl) . Bright sources have usually 
been assigned a spectroscopic redshift, with small error, while faint 
sources have been assigned a photometric redshift with larger error. 
Fortunately, even though these two different methods for redshift 
measurement have different error sizes, the infrared and X-ray lu- 
minosity distribution inferred both by spectroscopically and photo- 
metrically determined redshfits have the same median value. There- 
fore, we conclude the difference in redshift measurements does not 
induce any systematic bias. 

Figure[T]shows the redshift distribution in different infrared lu- 
minosity intervals, star-forming galaxies (SFGs; Lm. < IO^'^Lq) 
with median z ^ 0.168, luminous infrared galaxies (LIRGs; 
1O"L0 < LiR < IO^^I/q) with median z ~ 0.518, and ultra- 
luminous infrared galaxies (ULIRG; Lm ^ IO^^Lq) with median 
z ~ 1.268 (for more details on the computation of iiR please see 
section [331 l. 

^ http://www.cfht.hawaii.edu/~amouts/LEPHARE 
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Table 1. Samples properties and matching results 



Catalog 


Ai-ea 
in deg^ 


Flux limits 
in mjy 't) 
or lO^^^ergs cm~^ s~^ W 


Ns 




References f*^' 


Spitzer-COSMOS (S-COSMOS) 
70|xm galaxies catalog 


2.47 


8.5 (t) 


(5(t) 


1503 


1503 


Fraver et al. (2009) 
& Kartaltepe et al. (2010) 


XMM-COSMOS 
Point source catalog 


2.13 


7.2 (t) 
40 (t) 
97 (t) 


(Soft band) 
(Hard band) 
(Full band) 


1886 


108 


CaDpelluti et al. (2009) 
&Brusa et al. (2010) 


Chandra-COSMOS (C-COSMOS) 
Point source catalog 


0.50 (Deep) 
0.40 (Shallow) 


1.9 «) 
7.3 (« 
5.7 (t) 


(Soft band) 
(Hard band) 
(Full band) 


1761 


92 


Elvis et al. (2009) 


XMM- & Chandra-COSMOS 
concatenated catalog 








2821 


142 


This work 



'"'Effective area covered by the survey; Flux limit of the catalog; ''^'Total number of sources; ('''Number of sources cross-matched with the 70^m 
catalog; ("^'References describing the catalog; (^'Catalog concatenated from XMM- and Chandra-catalogs; sources matched between the two catalogs are 
not duplicated; the cross-match between this catalog and the 70|J.m catalog constitute our primary "70|xm/X-ray " catalog. 



Table 2. Confirmed spectroscopic redshift sample 



Telescope Instrument 



Reference 



ESO-VLT VIMOS 49 zCOSMOS 

Lilly et al. (2007) 

Magellan IMACS 26 Trump et al. (200T) 

Sloan SDSS 15 Abazajian et al. (2009) 

Keck II DEIMOS 11 Kartaltepe et al. (2010 ) 
MMT Hectospec 1 Prescott et al. (2006) 

("'Number of 70|J.m/X-ray sources with spectroscpic redshift. 



Redshift distribution 



SFGs (log{L,„)<ll) 
LIRGs (1 1<log(L„)<12) 
ULIRGs (log(L|f,)>12) 



■1 I 
I I 
L J 



3.3 Total infrared luminosity 

iKartaltepe etall j2010') estimated the total infrared luminosity 
(Z/ir) for all 70 Jim-selected galaxies, by fitting a SED with a x^- 
minimization method using Le PHARB^ and integrating the lumi- 
nosity over the 8-1000(i.m wavelength range. 

The current pu blic COSMOS 70iim-catalog from 
IKartaltepe et alj ( l201(]t) has been constructed carefully taking 
into account the possible AGN contamination: they only used 24, 
70, and 160 \im photometric data. However, this catalog suffers 
incompleteness in the 160|J.m band, and lacks any direct photomet- 
ric obs ervation in Ray leigh- Jeans part of the SED. IKartaltepe et al.l 
assessed that measurements of Lir without 160iim data are 
underestimated by 0.2 dex compared to measurements with 160n.m 
data. Despite these caveats, we estimate that this catalog is robust 
enough for our study, as the galaxies without 160nm photometry 
are equally spread between the samples with or without X-ray 
detected TOjim galaxies. 

We divided our 70|i.m/X-ray galaxy sample into three distinct 
classes according to their total infrared luminosity as star-forming 
galaxy (SF), luminous infrared galaxy (LIRG), and ultra-luminous 
infrared galaxy (ULIRG) with Lir < IO^Lq, 10"Lq ^ Lir < 
10^^ Lq, and Lir ^ 10^^ L©, respectively. The redshift distribu- 
tions of our three classes are represented in Figure |T] 



Figure 1. The redshift distribution of 70|J.m/X-ray galaxies. We divided 
our total sample into three distinct classes selected according to their total 
infrared luminosity : star-forming galaxies (SFGs; Lir < IO^^Lq; blue 
solid line) at a median redshift of z ~ 0.168, luminous infrared galaxies 
(LIRGs; 10" Lq ^ Lir < IO^^Lq; green dash line) at z ~ 0.518 and 
ultra-luminous infrared galaxies (ULIRG; Lir ^ 10^^ Lq; red dash-dot 
line) at z ~ 1.268. 



3.4 Estimating dust temperature of host galaxy 

In the local universe, the far-infrared luminosity is assumed to be 
driven by the emission from dust grains contained in the galaxy. 
Far-infrared photometry becomes a critical indicator to reveal the 
dust properties of the galaxy. Owing to the deep Spitzer observa- 
tions, we identified 463 (463 out of 1503; 31%) lQ\im sources, 
amongst which 52 (52 out of 142; 37%) of our 70|i.m/X-ray galax- 
ies, with 160iim counterparts. Most (89%) of 160[xm detected 
sources have a single matched 70iim source inside MIPS error cir- 
cle. While for the remaining sources with multiple 70n.m counter- 
parts, |K^alte2e|eral](|20l3) assigned the brighter one as counter- 
part of the 160[Am source. 

Using the lQ\im and 160iim fluxes, we adjusted a single tem- 
perature model represented by a modified blackbody radiation, as 
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following: 

F^ = v^B^{v,T) (1) 

with a fixed emissivity /3 ~ 1.5 l lClements et"Zll2010h 

We restrict the temperature range from 5K to 205K, typi- 
cal for dust temperature of submillimeter galaxies ( SMGs) and 
ULIRG jC hapman et al. |2005| : iKovacs et al] |2006| ; lYang et alj 
l2007l : ICasev et al.(20 09'). Models are convolved with Spitzer-MIPS 
filter response, and the best temperature is fitted using a minimum 
method: 



X 



E 



Fob 



■ nF, 



model, i -1 2 



(2) 



where i refers to 70nm and 160(i.m flux densities and n is the nor- 
malization factor to the 160nm flux density. 

To ensure an accurate estimation of the dust temperature in 
high-redshift galaxies, we exclude the 24|i.m fluxes from our SED 
fitting, mainly because of: 

(i) Contamination of the flux in the 5-30vim range by emissions 
from th e AGN dust torus or alternatively a p otential hot dust com- 
ponent (Brand et al.ll2006l : IPozzi et alJl201 2Y (ii) Presence of sev- 
eral Polycyclic Aromatic Hydrocarbon (PAH) transition lines (e.g. 
6.2, 7.7, and 11.3ixm emission lines) entering the 24^.m band at 
z~l, (iii) Presence of strong silicate absorption feature at 9.7|i.m, 
which is strongly correlated with the optical depth; in particular in 
the case of ULIRGs, presence o f large amount of dust grains causes 
the distinct absorption features dLevenson et al.ll2007h . 

3.5 Correction of X-ray properties 

With the help of combining both XMM and Chandra X-ray cata- 
logs, we extend our total sample of 70nm/X-ray galaxies to 142 
sources, of which 58 are listed in both X-ray catalogs. The count 
rates for these sources have been converted into 0.5 — 2keV (soft 
band) and 2 - lOkeV (hard band) flux density for both XMM 
and Chandra observat i ons by assumin g a power-low spectrum 
dCappelluti et alj |2009l : felvis et al.|[2009l) . For the 58 sources with 
overlapping detections between the two catalogs, we confirmed that 
both fluxes from XMM and Chandra are consistent, with a la error 
within 0.219 in log scale. 

Given the large redshift range covered by our sample (see Fig- 
ure [T), it is necessary to apply k-corrections to compute the X-ray 
rest-frame luminosities of our 70nm/X-ray galaxy sample. We as- 
sumed that the spectrum of our galaxies in X-ray follows a simple 
power law with photon index F = 1.7, and derived the ratio be- 
tween rest-frame and observe-frame flux integration, to compute 
the intrinsic rest-frame flux density. 

One of the main goals of our study is to explore the relation- 
ship between the galaxy host extinction and the central region ob- 
scuration. However, the neutral hydrogen column density identi- 
fied in X-ray is integrated on the line-of-sight, therefore estima- 
tions of the intrinsic colu mn density from X -ray spectrum fitting 
become model-dependent jAkvlas et alj2006l) . In order to simplify 
our analysis, we decided to use the hardness ratio (hereafter HR) as 
a good proxy of central region obscuration, define as: 

(3) 



H + S 



where H is X-ray hard band counts and S is X-ray soft band co unts. 

Our HR values are extracted from iBrusa et al.l (|2010|) for 
XMM observations and lElvis et ^ ( l2009h for Chandra observa- 
tions. To correct for the different observed bands between XMM 
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Figur e 2. Hardnes.s ratio extracted from th e XMM catalog jBrusa et al.l 
I2OIOI1 vs. value from the Chandra catalog l lElvis et aP |2009|) for our 58 
sources detected with both observatories. The different hard X-ray energy 
bands between the two observations, influence the hardness ratio measure- 
ments. We fit a simple line to estimate a conversion factor between the two 
different measurements. The fitted conversion factor is of 0.94 with small 
dispersion, consistent with the empirical expected value, predicted assum- 
ing a r = 1.7 X-ray power-law spectrum. 



and Chandra, we assume once again that our source X-ray spec- 
trum follows a power-law with F = 1.7. The final HR conversion 
factor from Chandra to XMM is ^ 0.94. To validate our empiri- 
cal conversion factor, we take advantage of our sources detected in 
both observatories. In Figure |2] we compared the HR from XMM 
with the HR from Chandra, for the overlapping detected sources. 
Just fitting a line between the HRs provides a similar conversion 
factor than our empirical estimation, proving the robustness of our 
conversion factor. Hereafter, all of our HR values are in XMM band 
system. 



4 RESULTS 

4.1 An AGN or star formation dominated sample? 

Although our sample is drawn from both 70\im and X-ray catalogs, 
it is possible that the X-ray luminosity is attributed to a combination 
of emissions from supemovae remnant, high mass X-ray binaries 
(HMXBs), and low luminosity AGNs, instead of being generated 
by a single powerful AGN. For the low X-ray luminosity galaxies, 
it is difficult to identify the origin of X-ray emission. For instance, 
Arp220 is considered as the typical ULIRG dominated by star for- 
mation. However, the detection of high ionization Fe K lines in its 
X-ray spectrum prevents us fro m completely ruling out the possi- 
bility of a low luminosity AGN jlwasawa et al J 20091) . 

Potentially, Ultraluminous X-ray sources (ULXs, X-ray point 
source; Lx ^ 10'^^ ergs~^) could explain the presence of X- 
ray emission in a star-forming galaxy. The presence of ULXs are 
mainly good indicators of recent star formation. However, deeper 
X-ray observations of Aip220 ruled out any ULX sources detec- 
tion, in contradiction to the current scenario of recent starburst. A 
plausible interpretation is that ULXs may be embed ded in dust, 
preve nting them from being detected in X-ray survey dSmith et al.l 
I2OI2I) . 

In this paper, we still define Arp220 as a star-forming ULIRG, 
as it is commonly done in the literature. Also, we are using Mrk23 1 
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as an archetype for AGN ULIRG, as indicated by its strong X-ray 
emission. 



4. 1. 1 AGN criterion from color-color selection 

Mid-in frared IRAC color-color selection described by IStern et al] 
( I2OO5L hereafter S05) is widely accepted as a robust method to se- 
lect the AGN candidates. Such color criteria are efficient to sep- 
arate AGN candidates from star-forming galaxies, due to the red 
[3.6^m] — [4.5iJ.m] color inferred by the combination of the power 
law continuum spectrum of the AGN and the relative weakness 
of the s tellar bump feature at 1.6iim usually emitted by the host 
galaxy. iLacv et al.l ( l2004h proposed another AGN diagnosis us- 
ing a different IRAC color-color selection from 805. However, 
IChoi et al. ( 201 1) used complete AGN samples, based on the differ- 
ent selections (e.g. BPT diagram, [O III]/H/3 ratio, high excitation 
line [Ne V], and broad line region feat ure), to show that the 805 
color-color selection is more robust than lLacv et al.l(l2004l) . 

In Figure[3] we show the IRAC [3.6(j.m]-[4.5nm] vs. [5.8iAm]- 
[S.Ojim] color-color diagram of our 70(i.m/X-ray galaxies. The 
deep pink line indicates the AGN dominating region defined by 
805. Using these criteria, we separated AGNs from our star- 
forming galaxies, LIRGs, and ULIRGs. Their estimated AGN frac- 
tion is 20% (4 out of 20; blue stars), 48% (29 out of 60; green 
square), and 92% (57 out of 62; red cross), respectively. There are 
63.4% (90 out of 142) of 70nm/X-ray galaxies classified as AGN 
from the 805 selection criteria. We remark also that ULIRGs are 
overlapping more with AGN region in the diagram than the locus 
of star-forming galaxies; this result is consistent with literature that 
demonstr ates an increasin g infrared luminosity as the AGN fraction 
increases jLee et al.l2010l) . 

Our galaxy sample is spanning a wide range of redshifts which 
may influences the color-color properties due to k-correction ef- 
fect. To better quantify this redshift effect, we placed the SED of 
Mrk23I and Arp220 in different redshifts, from z = to z = 3, 
and monitored the color changes in Figure|3] Templates for Mrk23 1 
and Arp220 are representative of AGN-dominated and star-forming 
dominated ULIRG, respectively. The Mrk231 colors never leave 
the 805 AGN region at any redshift because of its featureless 
power-law 8ED in mid-IR, while the Aip220 colors do not fol- 
low the 805 AGN criteria because that the different PAH emission 
lines entering into the IRAC photometry bands. The 805 criteria are 
therefore not extremely robust against Arp220-like star formation 
dominating ULIRG. 

4.1.2 Lx vs. Lib, relation 

The relation between hard X-ray luminosity and total infrared lu- 
minosity provides us with another method to determine whether 
star formation from host galaxy or AGN predominates the overall 
8ED. While color-color selection can be altered by the presence of 
specific lines (Donley et al. 2008), total infrared lumin osity is more 
robus t because it is derived from 8- lOOO^xm continuum dSmail et"al] 
l20Ilh . 

In Figure |4] we display the hard X-ray (2 — lOkeV) lumi- 
nosity against total (8 — lOOOvim) infrared luminosity for our 
70iJ.m/X-ray galaxies, according to their iiR, luminosity, subdi- 
vided into three sub-samples: star-forming galaxies, LIRGs, and 
ULIRGs (same symbols as in Figure[3). The X-ray luminosity ver- 
sus Infrared luminosity plot shows a tight correlation with Spear- 
man's p ~ 0.88. However, this strong correlation may be a mani- 
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Figure 3. [3.6]-[4.5] versus [5.8]-[8.0] color-color diagram of 70|J.m galax- 
ies with X-ray detection. The different Lir -selected galaxy sample follow 
the description of Figure[T] with blue stars, green squares and red cross in- 
dicating star-forming galaxies, LIRGs and ULIRGs respectively. The pink 
so lid continuous line represents the boundaries of the AGN region defined 
bv lStem et al Orange dash line and black solid line represent the 

evolution with redshift from 2 = to z = 3 of the colors of Arp220 and 
Mrk231 templates, respectively. The SOS criterion are not efficient in sepa- 
rating AGNs and Arp220-like star-forming dominated ULIRGs. 

festation of an observational bias. To confirm this, we need to es- 
timate the upper and lower limits in the Lx vs. Lir plot. In the 
case of the total infrared luminosity, we slice the whole 70|i.m cat- 
alog (Kartaltepe et al. 2010) into several redshift bins from z = 0.5 
to z = 3, and identify the maximum and minimum total infrared 
luminosity. For the hard X-ray luminosity, we simulate the lower 
X-ray luminosity in the same redshift bins using the Chandra flux 
limits, and measure the upper X-ray luminosity directly from the 
XMM catalog (Brusa et al. 2010). The results are displayed as gray 
regions in Figure |4] These observational limits and their evolution 
with redshift indeed explain the tight coiTelation between the in- 
frared and X-ray luminosities. A deeper (and less extended) X-ray 
survey from 4 Ms Chandra Deep Field South (CDFS) has demon- 
strated that some X-ray faint sources with bright total infrared lu- 
minosity are indeed m issing from our shallower X-ray observations 
dMullanev et alj20I2l) . 

By studying X-ray properties of star-forming galaxies, 
[Ranalli et al. ( 2003) demonstrated that hard X-ray luminosity is 
proportional to star formation rate (SFR). In fact, the most impor- 
tant contributors of X-ray emi ssion in such galaxies are the high 
mass X-ray binaries (HMXBs).^ Franceschini et al] ( l2003l) explored 
the X-ray properties of a sample of ULIRGs without AGN sig- 
nature. The conclusion of their work is that the X-ray luminosity 
and spectral shape of most ULIRGs are dominated by hot thermal 
plasma and X-ray binaries, originated in recent starburst region. 

On the other hand, because of the presence of the parsec-scale 
dust "torus" surrounding the accretion disk of the cen tral SMBHs, 
AGNs can also emit light in the infrared wavelength ( Gandhi et al.l 
[2009i) . Indeed, the ultraviolet and optical light emit by the cen- 
tral accretion disk i s absorbed by the dust and reemit in infrared. 
iGandhi et alj (12009') observed the core of nearby AGNs with un- 
precedented high spatial resolutions in both mid-IR and X-ray 
wavelength, demonstrated a strong correlation between 12.3vim 
and hard X-ray luminosities. These observations unambiguously 
demonstrated the existence of AGN drive mid-IR emission with- 
out contamination from the host galaxy. Based on the intrinsic 
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Figure 4. X-ray luminosity versus total infrared luminosity of 70|J.m se- 
lected galaxies with X-ray detection. Symbols are the same than in Figure|3] 
The X-ray luminosity error bars are derived from the estimated flux en'ors 
and the total infrared luminosity error bars are derived from the Icr prob- 
ability distribution of the SED fitting. Gray regions represent the areas 
covered by maximum and minimum luminosities in different redshifts bins 
drawn from the full X-ray selected and and 70|J.m selected samples. These 
observational limitations are inducing the apparent strong correlation be- 
tween Lx and Li^ observed in in COSMOS field. Black dot line and black 
solid line represent the Lx vs Ljj^ re lations for pure AGN and and star- 
formi ng galaxy samples respectively l lMuUanev et al.ll201 ll : iRanalli et alj 
l2003h . The luminosities of star forming submilUmeter galaxies (SMGs) and 
AGN SMGs, based on their X-ray s pectrum fitting, are represented by filled 
and open purple circles receptively l lLahd et al.lbOld) . Although 70|xm/X- 
ray galaxies suffer large contamination from star formation, their location in 
this Lx vs LiH diagram is consistent with this of AGN dominated SMGs. 



Figure 5. U pper panel : Luminosity-Temperature diagrain; black points 
and red points refer to all 70\im galaxies and 70|J.m/X-ray galaxies, with 
temperatures derived from a single black-body fit with fixed emissivity. 
Both total infrared luminosity and temperature error bars are estimated from 
the la probability distribution of the SED fitting. Temperatures depend 
strongly of redshift, given it was estimated from one photometric band only, 
which result in wanner temperatures at higher redshifts. Similar effects are 
observed with the SMG population, whose luminosity-temperature relation 
is shown as the cyan area, although they have typically lower dust tempera- 
ture for similar infrared luminosity (iChapman et a l. 2005). Lower panel : 
Average temperature for different bins of total infrared luminosity (average 
from the values of upper panel) in different total infrared luminosity inter- 
val; black points and red points refer to all 70|J.m galaxies and 70|J.m/X-ray 
galaxies. The presence of AGN does not enhance the dust temperature of 
70|a.m-selected galaxies. 



AGN/ quasar IR SED from iNetzer et al.l dlOOTl) . iMuUanev et"aLl 
Em l]) related the 12.3iim luminosity to the to tal infrared luminos- 
ity ( Lir) (their equation 5) . By combining the lGandhi et aljf2009h 
and iMuUanev et alj ( l201lh relations, we derived the unbiased Lx 
vs. Lir relation for AGNs, free of any contamination from the host 
galaxy. This "pure" AGN Lx vs. Lir, relation is represented by the 
dotted line in Figure|4] Such relation has been confirmed by swift- 
BAT X-ray selected AGN population ( lTuelleretal.l2010h . together 
with IRAS infrared measurements (BAT/IRAS AGN). Indeed, their 
60iJ.m luminosity is more likely powered by the AGN rather than 
star formation activity from the host galaxy because of linearly in- 
creasing correlation between Lir and Lx (Mullaney et al, .2012 ). 
For our COSMOS sample, the Lx vs. Lir distribution of 70|im/X- 
ray galaxies deviates from AGN relation and star forming relation 
by 0.5 — 1 dex and 1 — 2 dex, respectively. From the infrared and 
X-ray continuum perspective, AGNs approximately dominate the 
entire system. 

SubmilUmeter galaxies (SMGs) are a population of objects se- 
lected according to their detection in submillimeter wavelengths. 
This population is dominated by strongly star-forming galaxies at 
high reds hift, with this star form ation producing a large amount of 
cold dust jChapman et alj2()0^ . The concomitance in redshifts and 
infrared luminosities with ULIRGs indicates tha t SMGs may be an 
early stage of evolution of the merger scenario dGreve et al.ll2005l : 
[Biggs & Ivison 2008; Engel et al. 2010). Ultra-deep X-ray obser- 
vations for those distant star-formin g galaxies i n dicate that 20-30% 
of SMGs host an AGN (Alexander etai]|2005l) . lLaird et all l l2O10h 
have studied the X-ray spectral properties of SMGs to classify them 
as AGN or starburst. These objects are represented in Figure |4] 



with open and closed purple circles, respectively. Our 70iJ.m/X-ray 
galaxies appear to share the same Lx vs. Lir relation as the AGN- 
dominated SMGs. 

In summary, all indicators in the Lx vs. Lir distribution of 
our 70(im/X-ray galaxies converge toward a AGN dominated pop- 
ulation. 



4.2 Dust temperature of host galaxy 

According to the predictions from the merger scenario, galaxy in- 
teraction triggers massive star formation, generating large amount 
of dust, at the same time, it disturbs the gas distribution, funneling it 
toward the central black hole, feeding the a ccretion disk, and the re- 
fore generating a phase of nuclear activity dlmanishi et al.l2010l) . If 
the center of the galaxy harbors an AGN, dust in host galaxy could 
be influenced by the radiations produced by the nucleus. One way 
to trace radiative feedback from AGN, is to observe the galaxy in 
the far-infrared wavelengths, a s they are a good tr acer of the dust 
temperature in the host galaxy d Magdis et al.l201^ . 

In order to accurately determine the wavelength at the peak 
of blackbody radiation, and separate several cold dust compo- 
nents with different temperatures, information on the fluxes in the 
Rayleigh- Jeans tail of the blackbody, hence at longer wavelength, is 
essential. However, in a simple case, assuming a single blackbody 
profile with fixed emissivity /? — 1.5 in the Rayleigh- Jeans regime, 
JVoixm and Fl6o^m fluxes are sufficient enough to fit Wien regime 
of the blackbody. In this simplified case, we can derive the upper 
limit of dust temperature from the limited information (for details 
of the calculation, see section [33t . 
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Figure 6. Kolmogorov-Smiraov test: the hardness ratio cumulative prob- 
ability distribution of XMM, Chandra, and 70i^m/X-ray selected galaxy 
samples represented in black, red, and green color lines respectively. Dis- 
tribution of XMM and Chandra indicates these samples are drawn from the 
same population. However 70|xm/X-ray galaxies seems not drawn from the 
same population than purely X-ray selected samples. The distribution indi- 
cates an excess of HR at ^ -0.3, reveahng that the cold dust from the host 
starburst galaxy me be responsible of additional obscuration. 

As shown in Figure|5] the dust temperature of 70ii.m galaxies 
correlates strongly with the total infrared luminosity. This is actu- 
ally an observational bias because of a monochromatic wavelength 
for estimating the dust temperature. Indeed, at higher redshifts, the 
rest-frame flux corresponds to shorter wavelengths, and therefore 
probes warmer temperature. A similar incompleteness bias is act- 
ing on the estimation of dust temperature for SMGs, the selection 
at longer wavelengths (observed frame 850ixm, i. e. ~ SOOyim at 
z ^ 2) probed a lower dust temperature regime (IChapman et al.l 
l2005h . 

However, incomplete bias have little influence in our com- 
parative study between 70nm-selected with and without X-ray de- 
tection. In total, we have estimated the temperature for 463/1503 
(31%) of the 70tim-selected galaxies, of which 52/142 (37%) 
with X-ray detection. Such a large sample of 70 (xm galaxies 
with/without X-ray provides a good statistical ground to investigate 
any difference in dust temperature between AGN and non-AGN 
galaxies. 

In Figure |5] 70 (xm galaxies with/without X-ray detection are 
labeled as red/black filled circles, respectively, with the error bars 
estimated from the Ict probability distribution of fit. Although 
70ixm/X-ray galaxies display systematically slightly higher dust 
temperature compare to non-X-ray detected 70 [im galaxies (a shift 
of 1.85-4.25K for a temperature of 36.3-73.3K), the difference is 
not significant given the large uncertainties in the measurements 
(between 5.50K and 17.88K). Therefore, we do not observe any 
noticeable effects of the presence of AGN on the dust temperature 
of the host galaxy, in agreement wit h the recent observ ations based 
on Herschel far-infrared result from lElbaz et al.lfcOlOh . 



4.3 Hardness ratio of 70ii.in/X-ray galaxies 

To explain the diversity of AGNs (e.g., narrow-lines vs. broad- 
lines), a unified model, based on the variation of obscuration due 
to the orientation of the dust torus surrounding the SMBH accre- 
tion disk, has been proposed. In such model, the presence of an 



edge-on dust torus not only blocks the broad-line emissions, but 
also increases the absorpt ion in soft X-ray because of the higher 
hydrogen column density jAntonuccUl 19931) . 

In addition, for on-going galaxy mergers, the triggered star 
formation and higher dust production generate a higher neutral 
hydrogen column density as well, absorbing soft X-ray (e.g., 
lEsguei et al.l I2OI2I) . Indeed, using the Extended Chandra Deep 
Field-South (ECDF-S) dataset. lTreister et al.l ( l2009i) observed large 
obscuration of the soft X-ray emission in star-forming galaxies 
hosting an AGN at high redshifts, revealing a large amount of neu- 
tral hydrogen density on the line-of-sight. 

We like to investigate if there is a direct connection between 
the obscuration properties of the AGN and those of the host galaxy. 
As absorption affects the soft X-ray emissions more than the hard 
X-ray emissions, the HR is a good tracer of obscuration, with 
higher value of HR corresponding to larger absorption. To test 
whether our two samples. X-ray selected galaxies and 70ixm/X- 
ray galaxies, are consistent with being drawn from the same sam- 
ple (i.e. our null hypothes is), we apply Kolm ogorov-Smimov (K-S) 
test {Numerical Recipes, IPress et al .I1992I) . 

As an initial test, we examined the HR of XMM and Chandra 
X-ray selected sources, to test the robustness of our conversion be- 
tween XMM HR and Chandra HR (for details, see section [l4l l. The 
results of our K-S test for these populations are shown in Figure|6] 
with the HR cumulative probability distribution of XMM selected 
sample in black and Chandra selected sample in red, and the K-S 
parameters are summarized in Table|3] Our result is consistent with 
XMM and Chandra selected samples being drawn from the same 
population. Although two X-ray facilities have different sky cov- 
erage and depth, XMM HR distribution is identical to the Chandra 
HR distribution. 

Then, we apply the K-S test between the 70ii.m/X-ray and 
XMM samples. This time, the null hypothesis is rejected at a signif- 
icance of 99% confidence level, implying that 70ixm/X-ray galax- 
ies and XMM galaxies are drawn from different populations. The 
K-S test between 70|J.m/X-ray and Chandra selected samples also 
rejects the hypothesis. Figure |6] displays the cumulative probabil- 
ity distribution of HR from XMM (black solid line), Chandra (red 
solid line), and 70ixm/X-ray (green solid line) selected samples. 
The K-S parameters resulting from the three tests are summarized 
in Tablelll 

The tests indicate that 70ii.m/X-ray galaxies include more ob- 
scured AGNs than purely X-ray selected samples. This implies that, 
in the case of AGNs hosted by dust enshrouded galaxies, the pres- 
ence of an additional physical process not accounted in the AGN 
unified model scheme. We speculate that the excess of obscuration 
is not only attributed to absorption on the line-of-sight by contin- 
uum or clumpy dust torus around the SMBH, but also is associ- 
ated with an additional more diffuse dust component generated by 
strong star formation in the host galaxy. The starburst can be subse- 
quence of an event involving the whole system (e.g. merger) or gen- 
erated by the circumnuclear region (e.g. bar). That is to say, AGN 
is hosted by a galaxy, hence we can not neglect the influence of the 
host galaxy on AGN properties. Only X-ray observations with very 
high spatial resolution will enable us to identify the origin of the 
extra amount of neutral hydrogen. 
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Table 3. Results of the K-S test on the hardness ratio between different 
catalogs 



Number(red) 
5 10 15 



Catalogs 




^max 


Significant 
level 


Hypothesis 


XMM& 


Chandra 


0.0397 


0.466 


Not Reject 


XMM& 


70^m/X-ray 


0.1918 


0.0054 


Reject 


Chandra 


& TO^m/X-ray 


0.2105 


0.0017 


Reject 



5 DISSUSSION 

5.1 AGN influence on the host galaxy dust temperature? 



Our analysis in section [42] demonstrates that we have no strong 
evidence of any influence by AGN on the host galaxy dust temper- 
ature. Although o ur result is in agreement with recent result (e.g. 



Elbaz et it is in contradiction with the recent work by 



Raffertvetal] ( l201lh 



Using infrared color log( J24um/f7omT i ), to sample the 
blackbody temperature profile, iRaffertv et al.l fcOllh concluded 
that galaxies hosting an AGN present a high dust temperature. We 
reproduced their analysis for our sample, as shown in Figure |7] 
The 70ti.m/X-ray galaxies present a higher color index in average 
th an overall 70 [xm gala xies at all redshifts, a similar conclusion 
to IRaffertv et al] | |201 ih . Though the 7Q\x.mlX -ray galaxies show 
higher log(J24u. m/f70ii.m), they still agree with th e local defini- 
tion of cold dust dde Gt-ijp et al.lll985l : ISanders et al.lll9 88b). 

To explain the difference in color index, we looked at the evo- 



lution of the F' 



24(im/-FVoiim 



color as a function of the redshift 



for two typical objects, the star-forming ULIRG Aip220 and the 
AGN ULIRG Mrk2 31. We used the galaxy templates provided by 
IPoUetta et all j2007h . as shown in Figure |7] In the case of Arp220, 
the rest frame 9.7vim absorption line is shifted into the 24(im band 
at redshift z ~ 1.5, causing an apparent lower value of the color in- 
dex. Similarly, at 2 > 2, the 7.7 and 8.2 (xm PAH emission lines are 
entering the 24)i.m band, inducing a higher F24(j.m/J70Vi.m color 
ratio. In the case of Mrk231, the AGN driven power-law contin- 
uum produces a more stable color index. Our 70nm/X-ray samples 
follow the same color index of Mrk231-like galaxies across the red- 
shift range < z < 3, and the remaining parts of our samples fol- 
low more a Arp220-like profile. Therefore, the difference observed 
in F2A\im / F-jQ^ira color between the two samples has little to do 
with the average diffuse dust temperature. 

To have a robust measurement of the variation of dust tem- 
perature in the host galaxy, we need to explore the longer wave- 
lengths that have less perturbation from emission and absorption 
line features. For instance, measuring the longer wavelength color 
(e.g. JVoi^m/fieo iim) is a mo re accurate method to estimate the 
dust temperature (ICasevll2012l) . In order to simplify the temper- 
ature fitting, the models of temperature usually assume that j3 
is a constant value. However, it depends on which photometric 
bands are included; shorter wavelength photometry will induce 
the lower /3 v alue, and probe a systematically higher temperature 
jMagnelli et al.l2012 ). Our fitting procedure provides no strong ev- 
idence that the cold dust of host galaxy is heated by AGN (see 
section 14. 2t . Such result is consistent with t he latest work from 
the Herschel observations ( lElbazetal.ll20IOl) . Actually, the pho- 
tometry in the 250vi.m, 350nm, and SOOvim bands provides infor- 
mation on the Rayleigh- Jeans portion of the spectra, and therefore 
enables a more accurate estimation of the dust temperature of the 
host galaxy. Such precise measurements show a 2 — 3K difference 
between the dust temperature of a galaxy with AGN, and without 
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Figure 7. Distribution of (-F24|xm/^70|xm) color index for 70|xm selected 
sample and 70n.m/X-ray sample. Left panel : Distribution of the color 
against redshift for the 70\aa galaxies (black crosses) and the 70|xm/X-ray 
sample (red crosses). The black and orange lines indicate the evolution of 
Mrk23 1 and A rp220 flux ratio as a function of redshift, computed using the 
templates from jPoUetta et al.) j2007h . Right panel : Histogram of the color 
index for 70iJ.m-selected galaxies in black and 70n.m/X-ray in red, with the 
red histogram normalized to the peak of black sample for better display 
(bottom and top axis represent the number of 70|J.m-selected galaxies and 
70|xm/X-ray respectively). The 70|xm galaxies have a median color index of 
log(i<'24|xm/^70|xm) ~ —1.25, while for 70|xm/X-ray galaxies it is of ~ 
—0.99. Both satisfy the local cold dust definition of log(J24um/JVoum) 
^ -0.7 jde Grijp et al.lll985l:ISanders et al.ll988bl) . 



AGN i lElbaz et al.ll2oT ^. This implies that the far-infrared emis- 
sion of the galaxy is dominated by the dust blackbody emission 
heated by st ar formation, in agreem ent with previous studies. As 
an example, Isilverman et al.l ( |2009|) reached a similar conclusion 
from the optical zCOSMOS dataset. Using [O III]-corrected [O II] 
emission lines as star formation indicator, the correlation between 
SFR and X-ray luminosity, they demonstrated that there is ongo- 
ing star formation in galaxies hosting AGNs. With the exception of 
the warm dust compo nent (as shown by the f 24ti.m/J70ti.m color), 
neither this work nor ISilverman et alT(l2009h seem to display any 
truncation of the star formation by AGN radiative feedback. 

5.2 AGN obscuration from starburst? 

According to the major merger scenario of galaxy evolution, de- 
veloped to interpret a potential ULIRG-QSO connection, galax- 
ies undergo an obscuration phase after the merger occurred due 
to an enhancement of star formation and dust production. This ob- 
scured phase is thought to be ended by a consecutive AGN phase, 
during which feedb ack from the central SMBH repels the dust 
( iHopkins et alj2008l) . Giving an illustration, the NGC 6240 Seyfert 
galaxy is a well studied ULIRG, with high infrared luminosity in- 
dicating a starburst activity whi ch generates a large amount of dust 
and molecular gas Jlono et al.ll2007l) . A detailed X-ray spectrum 
from Chandra of NGC 6240 reveale d a very strong a bsorption from 
the neutral gas column density Nh JPtak et alj2003l) . Thanks to the 
higher spatial resolution, Chandra observations of NGC 6240 iden- 
tified the pres ence of a double AGN s system, hidden in the core 
of the galaxy jKomossa et 111120031) . These observations are con- 
sistent with the idea that galaxy was formed consequentially of a 
merger . 

To efficiently investigate the evolutionary stages of the AGN 
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concomitance with star foimation, our approach is to study 
the properties galax ies detected in both infrared and X-ray. 
iRaffertv et alj ( 1201 ih had a similar approach to ours, looking for 
the X-ray counterpart of 70vim-selected galaxies in the ECDF-S 
and EGS fields. Performing both a X-ray spectrum and a X-ray 
band ratio analysis, assuming a power-low spectrum with different 
column densities, they did not identify any excess of neutral hy- 
drogen in 70nm-selected X-ray sources compared to other X-ray 
selected AGNs. They concluded that the absence of obscured AGN 
among the star-forming galaxi es contradicts to the current AGN 
and galaxy co-ev olution model ( Hopkins et ai]|2008h . Actually, the 
mam reason whv lRaffertv et all J2OI ih " did not detect any additional 
obscured AGN from star-forming galaxies, is probably the statis- 
tical insufficiency of their sample. They only had accessed to 17 
of AGN detected in 70vim, with hydrogen column density Nh de- 
rived by X-ray spectral analysis dTozzi et al .I2OO6I) , comparing with 
whole AGN population. 

In contrast, our large sample (142 70[J.m/X-ray galaxies) en- 
ables us a proper statistical study. We therefore took a different ap- 
proach by applying a K-S test on the distribution of hardness ratio 
to investigate if 70iJ.m/X-ray galaxies and X-ray selected AGNs are 
drawn from the same sample (see section|43}. Trichas et al. ( 2009) 
attempted to conduct such study based on shallower 70 (xm observa- 
tions from the Spitzer Wide Area Infrared Extragalactic (SWIRE) 
survey, but did not have sufficient statistics to conclude if 70 (xm 
galaxies with X-ray detection were different from the global X-ray 
population. However, owing to the deeper 70(i.m observations in the 
COSMOS field, we could extract a larger sample, and our K-S test 
shows in less than 1% of probability for 70|J.m galaxies with X-ray 
detection to be drawn from the global X-ray population. Such re- 
sult is in agreement with the merger galaxy formation scenario: the 
neutral hydrogen obscuration in infrared-luminous AGN not only 
comes from the dust torus component around the AGN, but also 
from extreme star formation region in the host galaxy. 



AGN/starburst co-evolution model (iHopkins et al.ll2008l) . A K-S 
test on the modified HR between 70^Lm/X-ray galaxies and whole 
X-ray samples rejected the assumption of them to be originated 
from the same population within 99% confidence level, supporting 
the idea that the excess of X-ray absorption does not come from the 
AGN dusty torus obscuration, but probably from additional diffuse 
obscuration generated by the star formation in the galaxy. 

(iii) We estimated the dust temperature of our 70 (xm galaxy 
samples by fitting their far-infrared photometries with fixed emis- 
sivity models. Despite the presence of a warm dust component, the 
cold dust shows a similar temperature in host galaxy with and with- 
out AGN, indicating that the longer wavelengths are still dominated 
by star formation. This evidence conflicts with a scenario where ra- 
diative feedback from the AGN truncates the star formation in the 
host galaxy. 



ACKNOWLEDGEMENTS 

The work presented here is supported by the National Science 
Council of Taiwan under the grants NSC99-21 12-M-003-002- 
MY2, NSC99-2I19-M-003-005, and NSC 99-21 12-M-003-001- 
MY2. We would like to acknowledge the wonderful work made 
by Jeyhan S. Kartaltepe on the Spitzer-COSMOS 70iJ.m point 
source catalog on which our analyses are based. We would like 
to thank the following people for very useful discussions, sugges- 
tions, and comments: Lin-Wen Chen, Jeyhan S. Kartaltepe, Wei- 
Hao Wang, Albert Kong, Matthew A. Malkan, Nick Scoville, and 
David B. Sanders. This research used extensively the COSMOS 
survey archive data of the NASA/IPAC Infrared Science Archive, 
which is operated by the Jet Propulsion Laboratory, California In- 
stitute of Technology, under contract with the National Aeronautics 
and Space Administration. 



6 CONCLUSION 

We have investigated the properties of 142 galaxies both detected 
in X-ray and 70ixm in the COSMOS field. X-ray data are obtained 
from both XMM and Chandra point source catalogs, and 70 (xm 
photometry is drawn from Spitzer-MIPS 70 [xm point source cata- 
log. We classified our sample into three distinct subsamples accord- 
ing to their respective total infrared luminosity (Lir): star-forming 
galaxies (Lir < 10^^1/0), luminous infrared galaxies (LIRGs, 
10^^1/0 ^ I/iR < 1O^^L0), and ultra-luminous infrared galax- 
ies (ULIRGs, LiR ^ 10^^1/0), with median redshifts of z~ 0.168, 
0.518 and 1.268, respectively. The major conclusions for this study 
are as follows: 

(i) We applied two methods to determine which mechanism 
dominates the SED, star formation or AGN: 

i) Using Spitzer-IRAC colors, we have shown that the majority of 
our sample is dominated by AGN. Although the higher AGN frac- 
tion accompanies with higher total infrared luminosity (Lir), the 
presence of PAH emission lines related to star formation could dis- 
turb AGN identification. 

;/) Using the relation between hard X-ray (2-lOkeV) and total in- 
frared luminosity Lir (8-IOOOvim ), we found that our samples are 
dominated by AGN, with property comparable to AGN detected 
SMGs. 

(ii) We provided evidences for additional X-ray obscuration 
in X-ray detected 70ixm galaxies, in agreement with current 
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